Direct visualization of aging in colloidal glasses 
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Abstract. We use confocal microscopy to directly visualize the dynamics of aging 
colloidal glasses. We prepare a colloidal suspension at high density, a simple model 
system which shares many properties with other glasses, and initiate experiments by 
stirring the sample. We follow the motion of several thousand colloidal particles after 
the stirring and observe that their motion significantly slows as the sample ages. 
The aging is both spatially and temporally heterogeneous. Furthermore, while the 
characteristic relaxation time scale grows with the age of the sample, nontrivial particle 
motions continue to occur on all time scales. 
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1. Introduction 

Colloidal suspensions have long been used as model systems to study the glass transition 
in, H, H, As the concentration of a colloidal suspension is increased, the motion of 

the colloidal particles becomes increasingly slowed. The glass transition is considered 
to occur at the concentration at which particle motion ceases to be diffusive on long 
time scales . Recent work has studied the behavior of simple colloidal suspensions as 
this transition is approached |3, §]. In this paper we examine the behavior of glassy 
samples. Unlike supercooled colloidal fluids, the behavior of a colloidal glass depends 
on the history of the sample: particle motion becomes increasingly slowed as the sample 
ages [P, [l^, [Tl|, 0. We use confocal microscopy to study the three-dimensional 
motions of colloidal particles in an aging colloidal glass. We find that aging occurs 
due to slight rearrangements which are both spatially and temporally heterogeneous. 
Moreover, even after the sample has aged for long times, small but significant motions 
occur on short time scales. Previous work has seen similar evidence of spatial and 
temporal heterogeneities in aging systems [16], p!8|] . 
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2. Experimental Methods 

We use poly-(methylmethacrylate) (PMMA) particles of radius a = 1.18 fim and 
polydispersity ~ 5%, sterically stabilized by a thin layer of poly-12-hydroxystearic acid. 
The particles are dyed with rhodamine and suspended in a mixture of organic solvents 
(cyclohexylbromide and decalin) that closely matches both the density and index of 
refraction of the particles. We use a scanning laser confocal microscope to acquire 
images of a viewing volume of 63 fim x 58 fim x 12 fim at the rate of 3 images per 
minute. The viewing volume typically contains ~2400 particles. We focus at least 60 
fim away from the cover slip of the sample chamber to avoid wall effects. We identify 
particles with a horizontal accuracy of 0.03 /zm and a vertical accuracy of 0.05 fim and 
track them in three dimensions over the course of the experiment [|14|, |15|. 

The control parameter for the colloidal phase behavior is the sample volume fraction 
(f). While the rhodamine imparts a slight charge upon the particles, their phase behavior, 
4>freeze = 0.38 and (pmeit = 0.42, is similar to that of hard spheres {(pfreeze = 0.494 and 
4>meit = 0.545). We observe a glass transition at (pg ~ 0.58, in agreement with what is 
seen for hard spheres [0, We examined samples with volume fractions ranging (p ^ 
0.58 to (j) ^ 0.62. These samples form small crystals which nucleate at the coverslip, 
but do not form crystals within the bulk of the sample even after several weeks. 

In order to initialize the system, a small length of wire is inserted into each sample 
chamber. After placing a sample on the microscope stage, a handheld magnet is used to 
pull and rotate the wire through the sample for several minutes. The subsequent particle 
dynamics are reproducible after this stirring. Within a minute of ending the stirring, 
transient flows within the sample greatly diminish, and the particles move slowly enough 
to be identified and tracked. This defined our initial time t = for each sample, although 
the results below are not sensitive to variations of this choice. During the course of the 
experiments, no crystallization was observed within the viewing volume. Note that in 
many aging studies, the initial sample is prepared by a temperature quench (which 
would correspond to a rapid increase of the volume fraction in our experiments). The 
initial conditions in our experiments correspond to a shear-melted sample; the volume 
fraction remains constant. 

3. Results 

As a glass ages, we expect dynamical properties to change with time. To investigate 
this, we consider small temporal portions of a data set and calculate the mean square 
displacement (Ar^) (MSD) for each portion, where the angle brackets indicate an 
average over all particles and all initial times within each temporal portion. Figure 1 
displays the progression of the MSD curves of a colloidal glass (0 = 0.62) over a period 
of two hours, showing a marked change as the sample ages. The data is divided 
logarithmically into segments according to the sample age t^, the time waited since 
stirring, such that t^"*""^ = ct^, where c = 1.4, t° = 3.0 minutes and n = [1, 2, 11]. 
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Figure 1. Mean square displacement versus lag time for a colloidal glass {(j) « 0.62). 
The total length of the data set is 120 min. Data is divided logarithmically according 
to sample age. The geometric mean ages of the curves, from left to right, are tw — 3.0, 
4.2, 5.9, 8.2, 12, 16, 23, 32, 44, 62, 87 min. The dashed line has a slope of 1. 



The MSD for a given ty^ is calculated using data within the interval [c~^t^, ctw]- The 
MSD behavior at low and intermediate lag times is consistent with that of super-cooled 
liquids 0. For At < 30 s (not shown), the MSD grows roughly linearly with 
At, corresponding to the diffusion of particles inside "cages" formed by neighboring 
particles ||^, |^ . The MSD at intermediate lag times exhibits a plateau, as particles are 
locally confined by their neighbors. At larger lag times, the MSD curves each show a 
slight upturn, reminiscent of super-cooled hquids [l^. For super-cooled liquids, this 



upturn corresponds to cage rearrangements: the configuration of particles around a 
caged particle changes noticeably p. In our aging samples, we see similar, though 
smaller, motions, and the upturn in the MSD is less obvious than is seen in liquids. 
These motions, which change the local positions of particles, are likely related to cage 
rearrangements ||12|, [13| . Presumably these local rearrangements lead to the aging of the 
sample. However, unlike super-cooled liquids, the time scale for the upturn in the MSD 
depends strongly on the age of the sample, ty^, as is clear from Fig. 1 [jl3[. For larger 



values of t^, the plateau extends over a larger range of At. This behavior has been seen 



before in a variety of aging systems |10, 11 



We observe also that the MSD of a colloidal glass does not evolve uniformly. 
Consider the four rightmost curves in Fig. 1, corresponding to t^ = 32, 44, 62, and 
87 min. The first two curves (t^=32, 44 min) are virtually indistinguishable, showing 
that the dynamics have not changed within the imaging volume. However, the next two 
curves (t^=62, 87 min) are quite different, showing that the particles within the imaging 
volume have aged. This suggests that aging is in part temporally heterogeneous in our 
small observation volume. It is to be expected that if a larger volume was examined, 
such as in light scattering experiments, the dynamics would appear to evolve more 
smoothly with t^ [0, [ll[ . 



The temporal averaging used in calculating the MSD obscures some of the details 
of this temporal heterogeneity. However, it is unclear how to best study the dynamics 
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Figure 2. Various ensemble features as a function of sample age, tw The time lag 
is fixed (AT = 10 min) and the data are the same as shown in Fig. 1. (a) (Ax^) and 
7 = ■^/(A^^y/S (dashed). (The x component is used as the z component is noisier due 
to optical effects.) (b) Average number of particles in a cluster of mobile particles; see 
text for details. The dashed line indicates the average cluster size if mobile particles 
were randomly distributed. 





given that the time scales for such dynamics may vary dramatically with the sample age 
tw Furthermore, the mobile particles responsible for the upturn in the MSD are only 
moving slightly farther than the immobile particles; the latter still move locally due to 
their Brownian motion. To help distinguish the rearranging particles from those which 
only move within their cage, we average the trajectory of each particle over 1 min. The 
results that follow are not sensitive to the choice of this averaging time. 

To investigate the nature of the temporally heterogeneous aging, we choose a fixed 
lag time AT and study how the dynamical behavior changes with t^. We might expect 
that for tw -C AT, t^ ~ AT, and t^ S> AT, we should see strikingly different behavior. 
We choose AT = 10 min to allow us to resolve these three regimes with our data. 
The solid line in Figure 2(a) shows {Ax'^)at plotted as a function of sample age, 
here, the angle brackets do not indicate a time average, but only a particle average, 
unlike the MSD. As increases, (Ax^) generally decreases, which is not surprising 
given the behavior shown in Fig. 1: at large ty,, particle motion on a time scale of 
AT=10 min reflects the motion of particles confined by their neighbors. However, 
in addition to this overall relaxation of (Ax^), Fig. 2(a) shows a number of abrupt 
events, when {Ax^) jumps in value before relaxing further. Further confirmation of 
the temporally heterogeneous nature of this process is demonstrated by the dashed 
line in Fig. 2(a), which shows 7 = ^ (Aa;^)/3. If the dynamics were purely gaussian, 
i.e. diffusive, we would expect 7 = (Ax^). The consistently higher value of 7 and its 
prominent correlation to events in (Aa;^) indicate that the ensemble dynamics are highly 
nongaussian and events in {Ax"^) particularly so. 

It is intriguing that there is no unique behavior for t^ ~ AT; instead, aging "events" 
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Figure 3. (Ax^) as a function of sample age for various values of lag time, At, on 
linear axes (a) and log- log axes (b). The straight lines in (b) are fits to power law decay, 
(Ar^) ~ t^-'^, with a=0.23, 0.29, and 0.31 for At = 5, 10, and 20 min, respectively. 
The curves in (b) corresponding to At = 5 and 20 minutes have been vertically offset 
for clarity. The data are the same as in the previous figures. 



occur at many different values of t^j. To further confirm tliis, we plot (Ax^) versus sample 
age for several different values of AT in Fig. 3. The curves are clearly related; many 
events are simultaneously present in all curves. In Fig. 3(b) these curves are plotted 
on log-log axes, showing that the overall decrease in (Aa;^) is consistent with a power 
law decay, (Ax^) ~ tio~°. For different samples and different choices of AT, we find a 
ranges between 0.05 and 0.5. 

In super-cooled colloidal liquids it is known that the dynamics are spatially 
heterogeneous as well as temporally heterogeneous 0, §]. To investigate this for our 
aging samples, we plot the 3D locations of mobile particles in Fig. 4 for ty^ = 10, 55, 
and 95 min. The mobility of a particle is defined as its displacement during the interval 
tw to tw + AT (with At = 10 min). The particles shown are the 10% most mobile at 
those times, and are generally grouped into large clusters. Varying thresholds does not 
substantially change these pictures. Our results are similar to the spatial heterogeneities 
seen in simulations of aging [jl6[ . Signs of heterogeneities have also been seen in aging 
experiments which studied glycerol [|1^ . Previous microscopic studies of colloidal glasses 
did not show such large clusters |0; we find that the slight averaging of the particle 
trajectories discussed above (over a 1 min interval) is necessary to distinguish the 
rearrangements from the local Brownian motion of caged particles. Without averaging, 
our pictures look similar to those of Ref . ^ . 

It is striking that the mobile particles are clustered for all values of t^, notably at 
both tw = 10 min = AT and t^, = 95 min ^ AT. We investigate the t^ dependence 
of the spatial clustering of mobile particles by generating similar pictures for each t^ 
and measuring the average size of the clusters in each picture. (A cluster is defined 
as all connected mobile particles which are closer to each other than 3.2 fim, the first 
minimum of g{r), the pair correlation function.) Again, mobile particles are defined as 
the particles with displacements in the top 10% at each t^, using AT to characterize the 
displacements. The average cluster size shows large fluctuations, as seen in Fig. 2(b). If 
the 10% most mobile particles were randomly distributed in space, the average cluster 
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Figure 4. Locations of the 10% most mobile particles at three different ages t^. 
For each picture, mobility was determined by calculating displacements Ar over an 
interval [tyj,tu, + AT], with AT = 10 min. Left: = 10 min, and Ar > 0.43 /xm for 
the most mobile particles. Middle: = 55 min, Ar > 0.34 /xm. Right: = 95 min, 
Ar > 0.33 /Ltm. The data are the same as shown in previous figures, and the choices 
of tyj correspond to local maxima of 7 in Fig. 2(a). The particles are drawn to scale 
(2.36 fim diameter) and the box shown is the entire viewing volume (within a much 
larger sample chamber). 



size would be ~ 4.5 [dashed line in Fig. 2(b)], however, we find that the average cluster 
size is always larger than this value. The larger average cluster sizes correspond to 
locally large values of (Ax^) and 7, seen in Fig. 2(a). 

Our most surprising result is that overall, the cluster sizes do not increase 
dramatically with tw For the sample to age, particles must move, and may do so in 
small groups such as these clusters. For early ages t^j, it might be expected that particle 
rearrangements may involve very few particles, which can rearrange relatively quickly, 
thus causing the MSD to have a shorter plateau. For larger t^, if the rearrangements 
involved more particles, this would by necessity be more difficult and take longer to 
occur, thus explaining the slower dynamics reffected in the later MSD curves of Fig. 1. 
However, this is not what we see. As seen in Fig. 2(b), the average cluster size does not 
dramatically change as a function of ty,. This remains true even if we scale the choice of 
At with tyj (not shown). Thus, it is not yet clear what ultimately causes the dynamics 
to slow as our sample ages. While the rearrangements occur via groups such as those 
shown in Fig. 4, we have not found anything in their character which changes as the 
sample ages. 

While all the data shown in the figures correspond to one sample with (f) = 0.62, we 
have investigated five other samples in the range 0.58 < (p < 0.62. All data sets show 
similar behavior, and as yet we see no clear 4> dependence in our results. 

4. Conclusions 

Aging in our colloidal samples appears to be due to the most mobile particles in the 
sample, seen in Fig. 4 and reflected by the fluctuations in Fig. 2. In order for the 
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sample to age, particles must move to new positions. Our results show that such 
rearrangements are both spatially and temporally heterogeneous. In general, each aging 
event appears spatially uncorrelated with the previous. The surprising result is that 
significant fluctuations occur on time scales At ^ t^^,, despite the fact that for large t^, 
the mean square displacement shows only a plateau at such At (see Fig. 1). We find, 
in fact, that the plateau in the MSD is due to a temporal average over these infrequent 
rearrangement events and the more prevalent caged behavior. We have also attempted 
to investigate the behavior of the most immobile particles, but this has been difficult, as 
every particle still possesses nontrivial Brownian motion, and there is no unambiguous 
way to distinguish the most immobile particles. 

Thus, the changing character of the MSD (Fig. 1) and the related decay in (Ax^) 
for a fixed At (shown in Fig. 3) seem to be the most significant dynamical changes 
in aging colloidal glasses. In our future work, we plan to investigate the microscopic 
structural changes associated with these dynamical changes. 
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